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We examine how intense optical beams can penetrate deeply into highly absorbing media by a 
non-linear, photo-bleaching process. The role of stimulated recovery to the dye ground state can 
be important and is delineated. This analysis of non-linear absorption processes is applicable in 
general to situations where chromophores are irradiated, for instance in biology. We examine the 
implications for the bending of cantilevers made of heavily dye-loaded nematic photo-solids, that 
is nematic glasses and elastomers that have large mechanical reactions to light. In particular we 
describe the bending of cantilevers sufficiently absorbing that they would not bend if Beer's Law 
were applicable. We quantify the role of optically-generated heat in determining the mechanical 
response and conclude that in general it is minor in importance compared with optical effects. 



I. INTRODUCTION 

Nematic solids are principally elastomeric [l[ or glassy 
0. In each case they are composed of nematic polymers, 
but either loosely or tightly crosslinked. Their moduli 
are accordingly low or high and they are capable of ei- 
ther huge or modest extensions. They form unique solids 
in that, when their order is changed, they can change 
their dimensions considerably and these mechanical re- 
sponses can be steered in many ways. The order can 
either be redirected (director rotation under the influ- 
ence of electric, mechanical or optical fields) or reduced. 
In this paper, for simplicity we concentrate on chang- 
ing the magnitude of order, achievable by heating or 
cooling, or by the absorption of light into photo-active 
rods when they are a component of the nematic solid. 
The natural elongation of the oriented polymers mak- 
ing up the networks is modified by order change, with 
the result that macroscopic mechanical changes are in- 
duced. Reversible thermal strains of hundreds of per 
cent have been observed in elastomers [!, H|. Smaller 
changes arise in glasses, where subtle arrangements of 
director have been employed to achieve more complex 
mechanical responses |2j • Since light and heat have anal- 
ogous effects, one expects and indeed finds that optical 
absorption leads to analogous mechanical strains, both 
in photo-elastomers [f| 0] and in photo-glasses 0] ■ Such 
a mechanism of mechanical change offers the possibility 
of micro-opto-mechanical systems (MOMS) where ele- 
ments can be optically-induced to bend as elastomeric 
photo-swimmers [1] or as glassy cantilevers 0, H, fToj . 
The bend arises because the contraction is differential 
with depth since the optical beam is absorbed and is 
hence weaker with depth. How the beam intensity varies 
with depth because of non-linear absorption processes, 
and how this variation is translated into bending even in 
situations where absorption is so high that penetration 
would be negligible in the linear case, are the two ques- 
tions that concern this paper. A schematic, Fig.[TJ of the 
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FIG. 1: Radiation penetration with linear absorption length 
d giving light-induced bend. 



penetration and of the differential contraction leading to 
bend establishes the coordinate (x) of penetration, the 
thickness (w), and the curvature (1/-R) of the cantilever. 

Light disrupts nematic order when molecular rods 
which arc also chromophores, that is dye molecules, bend 
on absorbing a photon in their straight (trans ) ground 
state and make an indirect transition to their bent (cis ) 
excited state. The orientational order of tightly packed 
rods is thereby reduced by the cis fraction. We shall 
assume for simplicity that the mechanical contractile 
strain is proportional to the cis volume fraction. It is 
a vexed question as to whether heat released by these 
optically-induced transitions is the actual cause of the 
order change. This question has been addressed by ir- 
radiating polydomain nematic photo-glasses with polar- 
ized light. One finds [ll[ that the mechanical contraction 
occurs along the light polarization direction, a response 
that is therefore simply tuneable by rotating the light 
polarization. Domains closely aligned to the optical elec- 
tric field absorb, and thus also contract, preferentially 
over domains that are misaligned Moreover, since 

the domain size is small, heat diffusion is fast and do- 
mains neighboring an absorbing region quickly reach sim- 
ilar temperatures. If temperature were the mechanism 
for contraction, it would thus under these conditions of 
fast diffusion not distinguish a preferred direction and 
there would be no overall effect (just as when one simply 
heats a polydomain nematic solid). Evidently, optical ef- 
fects dominate. We discuss thermal effects at length in 
this paper and find their influence small, even across a 
whole monodomain cantilever, because of the diffusional 
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rates that arise. 

Since photons are absorbed by the dye components 
in the nematic, then light penetrating the face of a 
nematic photo-cantilever will be attenuated and hence 
the contractions generated diminish with depth. Cur- 
vature of the cantilever results, Fig. [TJ It is important 
in micro-opto-mechanical systems (MOMS) where ele- 
ments can be optically-induced to bend as elastomeric 
photo-swimmers H or as glassy cantilevers [jj • 

Weak beams decay exponentially with depth (Beer's 
law). The conversion of straight to bent (trans — > cis) 
forms of the dye molecules is also exponential in this 
limit. For a linear connection (valid for small strains) 
between cis population and contraction, in the Beer limit 
maximum cantilever curvature is predicted [13] for w/d ~ 
2.63, where w is the thickness of the cantilever and d 
is the exponential decay length: if w 3> d, only a thin 
skin of network contracts and its contractile stresses are 
insufficient to make the unstrained part of the cantilever 
below respond. Equally, if w <C d, then there is little 
variation of photo-strain through the thickness and the 
cantilever may contract but not differentially with depth 
and thus will bend little. The extent of bend, for a fixed 
w/d, was also predicted to be linear with intensity. 

Experimentally, cantilevers are commonly heavily dye- 
doped but still show appreciable mechanical effects [lj]. 
The above arguments suggest that this bending is unex- 
pected: high dye concentrations mean strong absorption 
and hence small d. In this thin-skin limit, d <C w, bend is 
expected to disappear. Evidently non-linear effects lead 
to deep penetration of light and thus lead to bending. We 
shall explain these related effects by making two assump- 
tions for simplicity to illustrate the underlying principles. 
We take the straight forward case of no photo-induced 
director rotation, that is for nematic glasses and con- 
strained systems (for instance where surface effects might 
be strong, where the director seems to be immobile under 
elongations imposed at an angle and also probably don't 
rotate during photo processes [15| . However, photoiso- 
merisation in glasses has additional subtleties; we derive 
the non-linear penetration of intense beams into glass as 
well. Secondly we assume that there is not appreciable 
reduction in the magnitude of the nematic order param- 
eter under illumination. This assumption is reasonable 
in glasses, but in liquid and elastomeric nematics, this 
assumption is too restrictive. Indeed a description of 
polydomain nematic elastomer response invokes both or- 
der reduction and rotation [121 ]. We return elsewhere to 
the role order reduction and rotation play in non-linear 
absorption and mechanics. 

An initial attack on this problem [161 ] invoked photo- 
bleaching, that is depletion of the trans form, letting light 
through to greater depths than would be expected by 
Beer's law. Here we derive these non-linear effects fully, 
paying attention to two additional influences that are po- 
tentially important. The optically-stimulated back reac- 
tion cis — > trans can alter the non-linear processes consid- 
erably and is of importance when the cis absorption line 



is not widely separated from that of the trans species. 
Secondly, thermal effects could be considerable, all the 
more so because of a possibly nearby nematic to isotropic 
transition where giant thermo-mechanical effects are 
known to occur, especially in elastomers. 

The remainder of this paper is organized as follows. In 
section [II] we derive equations which describe the attenu- 
ation of light passing through a region with photo-active 
chromophores. Our analysis of non-linear absorption is 
thus relevant to a wide range of situations where dye is 
irradiated, and is not limited to mechanics, which is our 
ultimate aim here. Our methods and results are similar 
to a penetrating analysis, experimental and theoretical, 
of Statman and Janossy [17j . In section [TV] we show how 
optical attenuation leads to cantilever bending and calcu- 
late the radius of curvature as a function of the incident 
flux of light. In section [V] we investigate the distribution 
of strain throughout the bent sample, and in particular 
we demonstrate that there can be several planes within 
the cantilever on which the net strain is zero. In section 
IVII we discuss the possible effects of temperature change 
owing to absorption of photons on the results presented 
thus far. We find that, for intense illumination, temper- 
ature distributions are symmetric about the cantilever 
mid-plane and hence do not contribute to bend, only con- 
traction. Finally in I VIII we present our conclusions. 

II. ATTENUATION 

We consider the situation of Fig. [2(b) of a long, slen- 
der cantilever of thickness w illuminated by light with 
incident flux Iq. The absolute number density of chro- 
mophores is p p h- At a time t after the onset of illu- 
mination and depth x within the cantilever, the frac- 
tion of these chromophores in the straight trans state is 
rit(x, t) and the fraction in the bent cis state is n c (x, t) 
I — nt(x,t). The magnitude of the Poynting flux at x 
and t is I(x, t). The dynamics of the trans fraction is de- 
termined by three processes, (i) an optically stimulated 
trans — ► cis reaction with rate Til(x,t)nt(x,t), (ii) an 
optically stimulated cis — > trans back-reaction with rate 
T2l(x, t)n c (x, i) and (iii) a spontaneous, thermally acti- 
vated, cis — » trans back-reaction with characteristic time 
t. r t and r c subsume absorption cross sections per chro- 
mophore and the quantum efficiencies $t c and $ ct of the 
stimulated trans- cis reaction and cis-trans back- reaction 
respectively, see [l?} for the rate equations in full with 
such factors explicitly given. We take the absorption 
cross sections to be independent of nematic order - as 
discussed above, changing nematic order is itself another 
source of non-linearity. Combining these three rates we 
obtain for the rate of change of the trans fraction: 

^ = -r t I(x,t)n t (x,t) + (^+TJ(x,t)j n c {x,t) (1) 

In this paper we confine ourselves to the steady-state, 
that is = 0. Setting this condition in eqn JTJ) and 
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taking out a factor of r gives the steady state trans and 
cis populations: 



n t (x) = 



i+r c ri 
i+(r,+r c )r/ 



( T ) - IVr/ 

cW - i+(r t +r c )r/ 



(2) 



where / is now I(x) simply, and will be determined below. 
We can identify two characteristic, material intensities, 
It = l/(r t r) and I c — l/(r c r). It is convenient to scale 
the flux by its incident value, thus l(x,t) — I(x,t)/Io- 
The reduced intensity is thus X = 1 at the entry sur- 
face x — 0, see Fig. Gib). We also define dimensionless 
quantities a = Io/It and f3 = Iq/I c - These are the in- 
cident flux reduced by fluxes it or I c characteristic of 
trans or cis molecules respectively. We refer to a and (3 
subsequently as reduced incident fluxes. Ignoring for the 
moment attenuation, a measures how much a beam in- 
tensity Io leads to trans conversion, eqn |T]), by compar- 
ing Jo to It , that is the ratio of the forward rate to the 
thermal backward rate, a = Io/It = -^0^/(1/7"). Like- 
wise, P is the ratio of the induced to the thermal back 
rates. 

In terms of the reduced incident intensities a and /?, 
the steady state trans and cis populations are given by: 



nJr\ - i+M, 



n c (x) 



ttl 

l + (a+/3)Z' 



(3) 



Here X is just X{x) since we have the equilibrium case. 
In the Eisenbach experiments [l8[ the average conversion 
was n c ~ 0.84. His measurements of attenuation sug- 
gested (3 ~ 0, and thus one can conclude from eqn ([3]) 
that his a ~ 5. Note that a and (3 are independent of 
chromophore concentration, but do depend on the choice 
of the light polarization [T3|. Experimentally it is easi- 
est to determine a for a system dilute in chromophores, 
where one can ignore the complications arising when at- 
tenuation is significant. These estimates for a are lower 
bounds on actual values; including the effects of attenu- 
ation through the cantilever will lead to higher values of 
a. In later work [j| one can deduce that a ~ 0.8. 

The divergence of the Poynting flux, at any point 
through the cantilever is equal to the amount of en- 
ergy taken out of the beam per unit volume per unit 
time. For simplicity, we ignore curvature leading to obliq- 
uity factors for the intensity of light falling on the sur- 
face. That is, we consider small deflections or diffuse 
light. Energy is taken out of the beam both by the opti- 
cally induced trans^cis reaction and by the stimulated 
cis—ftrans back-reaction, terms (i) and (ii) above. The 
divergence of the Poynting flux is thus related to the sum 
of the rates of these two processes. Thus: 



dl_ 

dx 



-^ t Y%I{x,t)n t {x,t) -j c TJ(x,t)n c (x,t) (4) 



where the constant 7 in each case subsumes the energy 
of an incident photon, fiw, the reciprocal of the quantum 
efficiency <f> for the relevant transition, and the absolute 
number density of chromophores, p p h (which could differ 
between trans and cis forms since bent molecules pack 



less efficiently). The appearance of <fr as an inverse is 
required since for each successful transition in the rate 
TJni (i = t,c) there will be unsuccessful absorptions that 
do not contribute to dnt/dt in (fTJ) , but nevertheless still 
deplete the optical beam and contribute to dl /dx in ([4]). 

Eqns and ^ are a pair of coupled, non-linear, first 
order partial differential equations for I{x, i) and nt(x, t). 
Solving these equations subject to the boundary condi- 
tions l(0,t) — Iq and nt(cc,0) = 1 is in general complex, 
though analytically possible in some limits. We return 
to this problem elsewhere and here take the time- 
independent, equilibrium state. Thus, using n c = 1 — ri t , 
dividing through by the incident intensity Iq, and let- 
ting dt = f/(7tr t ) and d c = f/(7 c r c ) denote the charac- 
teristic lengths for optical attenuation by trans and cis 
chromophores respectively, one obtains [HI [l7j 
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(5) 



In terms of the parameters a and (3, the ratio of the 
trans and cis lengths is dt/d c = (7c / 7t) {(3 / a) . The ra- 
tio r\ — 7c/7t is the ratio of the quantum efficiencies 
for the trans — > cis and cis — > trans reactions, that is 
i] = &tc/$ct- We shall take 77 = I in the numerical illus- 
trations in this paper. We have ignored any attenuation 
by the host material; one could include such effects by 
adding a simple Lambert-Beer term —I(x,t)/dh to the 
RHS of eqn ^ or eqn ([5]). See section [III Al for an analy- 
sis. 

Inserting the steady-state expression for n t from 
eqn ([3]) into eqn ([5]), and then integrating w.r.t. x, sub- 
ject to X(x = 0) = 1, we obtain (see also pL7|| ) : 



In I 



VP 



ft' 



In 



1 + /3'X 
l + (3' 



x 

"dt 



(6) 



where j3' = /3(1 + 77). In general this expression is very 
different from Beer's Law, X(x) = exp(— x/d). 

Deviations from Beer's Law come about because at 
high intensities the cis population increases (bleaching) 
and is generally less absorbing than the trans species. 
Optical penetration is then more effective and is of great 
significance for photo-mechanics. To most simply see how 
non-linearities (bleaching) manifest themselves, consider 
the limit j3 — > [l|| that arises when stimulated cis back- 
conversion is weak (for instance in the work of Eisen- 
bach). Under those circumstances r c ~ in eqn ([1]), and 
then takes the form: 



dl 
dx 



X(x). 



(7) 



A non-Beer form arises because n t = 1/(1 + aX) itself 
depends on X, eqn ([3]) . Integration gives 



In I + a{X — I ) = —x/dt 



(8) 



also a (3 — > limit of ([H}. A formal solution of which is 
X(x) = — Wi,(ae a ~^), where Wl(x) is the Lambert-W 
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function [2fJ. The non-exponential decay persists until 
around aX < 1 , whereupon n t becomes independent of X 
and reverts to simple exponential form. 

The limiting cases of absorption are important. 

(i) (a + (3)X <C 1. Now n t ~ 1 and n c « which ren- 
ders (O trivially of the Beer form. The limit obtains 
when a = Io/h and = Iq/I c are both <C 1 (since X is 
bounded by 1), that is, the incident beam is weak com- 
pared with both material fluxes It and 7 C . It also obtains 
when a and are not small, but when X <C l/(a + /3), 
that is when the beam has diminished (albeit linearly 
rather than exponentially, see the sketch below and also 
the high intensity traces of Fig. [2]) to the point that 
nt = const « 1 and Beer behavior is finally recovered. 
From ([8]) one sees in fact X(x) ~ exp[— (x — dta)/dt] for 
x > d t a, a shifted Beer form. Exponential decay is the ul- 
timate fate of all optical beams provided that cantilevers 
are thick enough to get this diminution of intensity. 

(ii) The high flux limit al 3> 1 (with ~ 0) is where 
the forward reaction dominates over thermal back re- 
action. Photo-equilibrium is highly biased away from 
trans, that is in eqn (|3"|) ru . l/(aX). Eqn (J7]) reduces 
to X' ~ — l/(ae?t), whence X ~ 1 — x/(ad t ). This linear 
penetration for x < adt is evident in Fig. [2] at the higher 
a values and is at the heart of why even highly absorb- 
ing systems can be responsive. When x > adt, decay is 
again exponential, see above. 

(iii) The high flux limit, 01 > 1, that is IT C > 1/t, is 
where the stimulated back-reaction dominates over the 
thermal back-reaction. In that case nt —> 0/(a + 0) and 
n c — > a/ (a + 0) are again constants and again ^ takes 
a Beer form in the equilibrium limit: 

dI/dx = -X0(l + rj)/(d t {a + 0)). (9) 

The decay is exponential, X = exp(— x/d e s) with an ef- 
fective decay length d e ff = d t (a + 0)/{0{l + r/)). 
Thus, except for /3 < 1, profiles start in a Beer-manner, 
have an intermediate non-exponential behavior if we have 
a > 1 with /3 C a, and then conclude with another Beer 
decay. The intermediate regime, 01 *~ 1, is where the 
stimulated back reaction rate is comparable to the ther- 
mal rate. 

(iv) When the decay lengths accidentally coincide, d t = 
d c = d, that is 7tT t = 7 c r c , one can easily see in ei- 
ther or in ([5]) that a Beer form pertains at all in- 
tensities or depths into the photoisomerizing medium: 
X = cxp(— x/d). 

Statman and Janossy fl7| investigated photo- 
isomerization of solutions of the commercially available 
azodye Disperse Orange (D03). They obtained a ratio 
of a. 10 ~ 5 while the accessible range of a is up-to ~ 80 
(corresponding to an incident flux of 15mW/mm 2 ). The 
spectral bands for trans — > cis and cis — > trans overlap 
quite strongly for D03; one might thus expect much 
smaller ratios of 0/a to be accessible when using dyes 
with more separated absorption bands. Indeed, Eisen- 
bach's attenuation study showed that, for his systems, 
d c ^> d t and thus that a 3> 0, possibly a ~ 100 x 0. 



We thus show results initially for high ratios a/0 for a 
range of incident intensities a and then look at smaller 
ratios where the non-linear region is not so pronounced. 

As can be seen from eqn ([3]) , the cis population is al- 
ways reduced when is finite. Thus for non-zero we 
require a larger value of a to achieve a particular value 
of n c . Similarly attenuation will lead to reduced inten- 
sities lower than unity in eqn ([3]), again requiring larger 
values of a to achieve the same cis concentration. Thus 
estimates of a from absorption are lower bounds on true 
values. 



A. CASE 1 - a/0 = oo 

We recall the case in which the illuminating light 
doesn't excite the cis — > trans back-reaction at all, i.e. 
= 0. Thus the intensity is described by eqn ([5]). Plots 
are given in [lj|, but Fig. [2] for the case a/0 — 50 of rel- 
atively small also shows how for low a the decay from 
the surface intensity is exponential, but that penetration 
is much deeper for higher a, being initially a linear decay 
until finally decaying as an off-set exponential beyond the 
characteristic depth dt, that is for x > d t a. 

Such non-exponential behavior suggests caution when 
trying to establish an extinction length from the atten- 
uation of a light beam on traversing a cantilever. For 
instance [1J] an attenuation of 99% on traversing a can- 
tilever of thickness 1/im, were Beer's Law being fol- 
lowed, would result from an extinction length of deccr = 
ljUm/(2 ln(10)) ~ 0.22/im. But if a is not small, much 
more light penetrates to x — w (Fig. [2] is a guide) . The 
cfaeer value derived is a gross over-estimate of dt ■ Solving 
for dt from eqn ([5]) for a given attenuation X(w) on reach- 
ing the back face at x = w yields for 99% attenuation: 
d ~ w/(a + 4.6) = l/xm/(a + 4.6) -> 0.04^m for a = 20. 
The true dt associated with a possible exponential decay 
may thus be much lower than the value e?Bccr estimated as 
above, an indication that light has penetrated much fur- 
ther into the sample than would be expected for a simple 
exponential profile. Quantitative measurements of light 
attenuation varying with thickness or varying with in- 
cident intensity would resolve this ambiguity about dt 
and also allow a determination of a. Attenuation at one 
thickness can only give an upper bound on dt- 

The reasons for departures from Beer's law for the in- 
tensity can be seen by returning the solution X(x) to 
eqn §3§ to obtain the spatial variation of the cis con- 
centration. Fig. [3] for the case a/0 = 50 is a guide to 
the = plots (see 16]) and displays exponential decay 
in n c (x) following X(x) for low intensity (a = 0.1, 0.5). 
High incident intensity not only lifts n c (x = 0) at the 
surface, but also flattens the decay with depth - high 
n c means low n t = 1 — n c and hence fewer trans dye 
molecules in a state to deplete the incoming beam (a 
photo-bleached state). With photo-bleached surface lay- 
ers, radiation penetrates well beyond x <~ dt, and equally, 
contraction extends deep into the bulk, certainly beyond 
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the Beer penetration depth dt- 

For a = 0.5 a point of inflection first appears at the 
surface, and moves inwards with increasing incident in- 
tensity a. The cis fraction at the inflection is always 
ri c = 1/3 in this model. In the general case where (3^0 
we find that the value of a at which the point of inflection 
first appears at the front surface is a complicated function 
of the constant ratio 0/a, and the value of the cis fraction 
at the inflection point is no longer 1/3. From eqn Q, the 
surface cis concentration is n c (0) = a/(l + a), since there 
1=1, and rises to saturation, n c = 1, as intensity in- 
creases. The precise form of the cantilever bend depends 
critically on the shape of these n c {x) curves. In particu- 
lar the development of a point of inflection allows three 
intercepts of the straight, geometric strain curves, and 
thus three neutral surfaces, we will later see. Since n c (x) 
changes shape with increasing intensity, we will find an 
elastic response that is highly non-linear with intensity. 
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FIG. 2: The decay in reduced light intensity with reduced 
depth for various reduced incident intensities a, with a weakly 
stimulated back reaction such that a//3 = 50. The plots are 
similar to those for /3 = 0, with small differences we discuss 
when considering heat generation. 



B. CASE 2 - a/p = 50 

We now consider the case with a weakly stimulated 
back-reaction, (3 = a/50, and henceforth take r\ = 1. 
Fig [3 shows plots of the reduced intensity as a function 
of x/dx- The reduced intensity curves are largely iden- 
tical to the (3 = case. Once again, for small values of 
a (=0.1, 0.5) the decay is essentially exponential, with 
a characteristic length given by d t . Increasing a leads 
to deeper penetration, with the initial decay being es- 
sentially linear. However close inspection of the a = 10 
curve reveals a slight upwards curvature, a result of the 
higher order corrections in (3 that take us from ([8]) to 
([5]) . Eventually T{x) becomes exponential at penetration 
depths x r*> d t a significantly greater than dt- 

As we show above, for (3 > 1, that is here a > 50, 
the initial behavior should (briefly) revert to being expo- 
nential before attaining the linear decay associated with 
non-Beer. This point is easier to display below when we 
consider smaller a/ (3 ratios. 

The cis profiles as a function of depth are shown in 
Fig [31 For the values of a plotted, the curves are essen- 
tially identical to those in the (3 = case. High incident 
intensities result in larger cis fractions near the surface, 
and a flatter decay as before. The point of inflection oc- 
curs for a cis fraction less than the 1 /3 found in the (3 — 
case, but it is still of importance that inflections exist for 
the character and number of neutral surfaces we explore 
later. 



C. CASE3-a//3 = 5 

Increasing the stimulated back reaction further we take 
a/ (3 = 5. The curves for the reduced intensity as a func- 
tion of depth, see FigfJl now differ somewhat from those 
in the small (3 limit, e.g. as in Fig [21 For the smaller 
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FIG. 3: The decay in cis number fraction n c (x) with re- 
duced depth for various reduced incident intensities a, with 
a /fl = 50. Increasing a extends the conversion to cis to 
greater depths because of photobleaching of the surface lay- 
ers. Reduced geometric bend strain for a = 5, 10 are 
shown as straight dotted lines. For cantilevers of thickness 
w — w* — 9.313d and w = 12. 5d respectively, there are three 
and two intersections (neutral surfaces) with the photo-strain 
curves. 



values of a the curves remain exponential with a charac- 
teristic length dt- Increasing a leads to some increased 
penetration, without showing the long linear decay in 
reduced intensity seen in Fig [2] For a = 10 we have 
(32 = 2 at most, a value evidently insufficient to sat- 
isfy limit (iii). We do not have a finite initial region for 
small x where the decay is exponential with d c Q given 
in and below eqn @. For these values of a and (3 one 
would have d e ff = 3d t - The dashed line shows the in- 
finite a limit, and corresponds throughout its range to 
eqn ([5]), i.e. an exponential with characteristic length 
d cS = dt(a + (3)/(f3{l + 1])) = 3d t . Note that the initial 
of the a — 10 curve is close to that of the a = oo curve. 
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FIG. 4: The decay in reduced light intensity with reduced 
depth for various reduced incident intensities a with a//3 = 5. 
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FIG. 5: The decay in cis number fraction n c (x) with reduced 
depth for various reduced incident intensities a with a/j3 — 5. 



III. OTHER ABSORPTION PROCESSES 

At least two other influences are important in non- 
linear absorption. How does one deal with the absorption 
of the host material that is not a chromophorc. In prin- 
ciple this is the simplest form of absorption and for weak 
beams where both processes are Beer-like the effect can 
simply be divided out. We examine how this procedure 
translates to the non-linear case. Another, more com- 
plex effect is that of the host medium when it is a solid. 
The isomerisation processes can be biased by mechanical 
effects, as is well-known from linear experiments. These 
influence the non-linear absorption for intense beams. 



A. Beer-Lambert host absorption 

It will always be the case that the host for the chro- 
mophores will also absorb light. For simplicity ignoring 
the backreaction, j3 = 0, the absorption equation ([5]) be- 



dl 
dx 



dt 



1 

dh 



I(x). 



(10) 



The absorption length dt of course emerges naturally 
by measuring the absorption of a sample without dye, 
— ln(X lu ) = w/dh. Using n t = 1/(1 + aT) 1 and denoting 
the ratio of absorption lengths tp = dt/d\ u integration 
yields: 



1 f 1 + ^(1 + al w ) 
^ 1 + ^(1 + a) 



= -w/d e & (11) 



where the effective absorption length is l/d e a — X/dt + 
1/dh- The relation (jTTJ) for T w has close similarities to 
that for the profile with a back reaction, eqn ©. The 
limit ip — > of eqn (JTTJ) of no host absorption leads to 
the previous Lambert-W expression given below eqn |3|) . 
In the conclusions we discuss how host absorption can be 
allowed for in measurements. 



B. Absorption of intense beams in solid hosts 

The role of the host has been much studied in dynam- 
ical studies of absorption. When the host, typically a 
polymer or a network, is below the glass temperature 
then it is observed that there is more than one charac- 
teristic decay time for the cis population of guest chro- 
mophores. Eisenbach [2l[ clearly observed two times, the 
relative weight of the processes being dependent upon 
temperature. The slower process was analogous to that 
observed in a liquid host and becomes the only pro- 
cess present at elevated temperatures. The faster pro- 
cess was speculated to arise from the strained state of 
a cis isomer; a fraction of the molecules when in a bent 
state are in conflict with the solid matrix around them 
and thereby decay more quickly to the trans form. Other 
authors [22| suggest that there is a spectrum of relax- 
ation times corresponding to a range of environments 
that cis chromophores find themselves in. Eisenbach also 
found that, when above its glass temperature, rubber 
provides a liquid-like environment for its chromophore 
guests. There was only one relaxation time and this was 
comparable to that found when chromophores were dis- 
solved in liquids, a fact of considerable importance since 
elastomers are used as photo-mechanical media. 

Ignoring for simplicity photo-induced back reactions 
(as shown to be the case for instance in the Eisenbach 
study), the dynamical equation (UJ is modified in that in 
principle the forward rate T t could become a function its 
environment (denoted by p) and the thermal back reac- 
tion rate 1/r is certainly a function of p. The result is 
that the characteristic reduced intensity is now a func- 
tion of environment and we denote it by a p — T^t p Iq. 
The equilibrium trans number fraction for chromophores 
in the environment p is thus n\ = 1/(1 + a p I), see 
the original equation ([3]). Let the probabilities of find- 
ing a chromophore in the p environment be a p . Then 
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the total trans volume fraction of chromophores is n t = 
^2 p a> p /(l + a p X), and it is this n t that must appear in 



(12) 
(13) 



eqn (J): 




1 dX 

X dx 


1 ^ a p 

d t ^ 1 + a v X 
v 1 


w 


r x dx 


~d~t ' 




w 

~d~t = 


1 / 

— 1 o,\ + a 2 — a - 
a V 



l-g \ 
l+a 2 l J 



+ ln(J) - (1 -X) 



a\a 2 
a 

aia 2 



In 



1 + aX 
1 + a 



(14) 



The first equation is general, the second two have been 
specialised to the Eisenbach case of two decay rates and 
thus Oi\ and a 2 with weights a and 1— a respectively. Here 
a = aa\ + (1 — a)a 2 ; note that this is not an averaged 
a p . Again, this result has similarities of form with the 
(3^0 form © and the host absorption result (fTTj) . and 
is derived by the same kinds of integrations. 



IV. OPTICALLY-INDUCED CURVATURE 

Fig. [1] shows a cantilever with radius of curvature R. 
The geometric strain from bending is x/R + K, where R 
is the radius of curvature adopted by the cantilever and 
if is a mean strain, both to be determined for a given 
thickness w and illumination. Illumination changes the 
natural length of the sample, the actual strain with re- 
spect to this new natural length is x/R + K — e p which, 
if we further reduce x/R and K by the dimensionless con- 
stant —A connecting photo-strain e p and the cis concen- 
tration, we obtain x/R + K + n c (x) for the effective re- 
duced strain. A linear relation between cis concentration 
and strain is probably valid for nematic glasses, but for 
elastomers it is possible to reach the isotropic state by il- 
lumination at temperatures close enough to the nematic- 
isotropic transition and the relation is no longer linear, 
but can be mapped on to the observed variation of strain 
with temperature The mechanical stress a is related 
linearly to the strain via the Young's modulus E. Since 
there are no external forces nor external torques, me- 
chanical equilibrium requires vanishing total force and 
moment across a section, thus: 



a(x)dx = E j (^ + K + n c (x)jdx = Q, 
xa(x)dx =E j x ^ + K + n c (x)^jdx = 0.(15) 



When the modulus is constant it drops out, but must 
generally be retained (for some photo-glasses E is known 
to vary with illumination [HI]). Performing these inte- 



grations we have: 

,,2 



W 

2R +KW 
w 3 Kw 2 

m + 2 



n c (x)dx 

i 

xn c (x)dx. 



(16) 
(17) 



Simplifying between these two expression we obtain for 
the radius of curvature: 



i_ _ 12 r /w 



x ) n c (x)dx 



(18) 



Eqn (O can be rearranged to give an expression for n c (x), 
recalling d t /d c = rj/3/a, r\ = 7 c /7t, and n t = 1 - n c : 



n c (x) 



1 dX 



l- v l i- v lXdx 



(19) 



Inserting this expression into eqn (fTS)) and changing in- 
tegration variables J^ 1 dx — > Jx™ = i dX we have: 

3 -, „x 



| = 12 [* 
R \w 



1 



l-r?£ 

' a 



W 

2dt 



x \ dX 



(20) 



Substituting for x/dt and w/dt from eqn ([6]) and inte- 
grating, we obtain (with (3 1 = /3(1 + 77)): 



R 



1 2it r 

-——3 Li 2 {-^)-Li 2 {-(3'X w ) 
/3'{w/d t ) L 

-iln(X w )ln[(l + /?X)(l + /3')] 



(21) 



where Li 2 (x) = f dt 1 ^ 1 '^ — J^kLiW ^ s ^ ne diloga- 
rithm [23]. The limit (3 — > within this expression re- 
covers our earlier expression for the curvature [161 ]: 



dt 
R 



12ad t 3 



(1-X U 



w 

2dt 



(22) 



At low incident light intensity, a — > 0, analysis [1J] 
for exponential decay gave maximal reduced curvature 
w/aR for w/d ~ 2.63. In this limit l/R ~ a ~ 7 , 
hence the division by a to obtain results universal for all 
(low) intensities of incident light. As intensity increases, 
the maximum in w/aR moves to larger w/d because the 
radiation penetrates more deeply. 

The non-linear regime, at fixed Beer's Law penetration 
d t , is best revealed by reducing R by d t instead of by 
w, and by not reducing l/R by a. Fig. [6] plots dt/R 
against w/dt to reveal maxima in d/R at greater w as 
intensity a and thus penetration increases. At a given 
w/dt, one sees curvature increase initially with intensity, 
the Beer limit, and then reduce as penetration increases 
and gradients of strain are reduced. Thus appreciable 
curvature arises experimentally even in cantilevers thick 
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FIG. 6: Curvature reduced by 1/dt against reduced cantilever 
thickness w/dt for various incident reduced light intensities a, 
with a//3 — 50 (upper) and = 5 (lower). For high intensities, 
bleached surface layers let light penetrate more deeply and 
hence a significant fraction of the cantilever has its natural 
length contracted. Bend occurs even for cantilevers much 
thicker than the linear penetration depth d t . Maxima in cur- 
vature occur at smaller thicknesses than in the = case 
where back reaction is purely thermal, and are seen at ever 
smaller thickness as f3 increases. 



in the sense w ^> d t : in [11] it appears that curvature is 
induced even though the cantilevers involved (with w = 
10/im) are apparently at least 100 times thicker than their 
extinction length! (See also [14] •) The curvatures against 
thickness for various intensities a in the case of a/ (3 — 50 
studied above are practically identical to those in the — 
limit of [l6| . As the back-reaction rate increases relative 
to the forward rate, a/f3 = 5, the penetration is less 
and the curvature maxima move to noticeably smaller 
reduced thicknesses w/dt, see lower pane. 

In the limit of intense illumination, that is a > 1, the 
intensity profile simplified to X ~ 1 — x/(ad t ) for pene- 
tration x < adt and hence intensity such that aX > 1. 
Under these circumstances n c ~ 1 — l/(aX) 2 and, impor- 
tantly for the curvature, n' c ~ — l/[d t (aX) 2 ]. If the thick- 
ness w -C adt j then the largely linear profile of X(x) ob- 
tains through most of the thickness of the cantilever and 
also I ~ 1. One can then approximate n' c ~ — \/{dta 2 ). 
While the strain follows n c (x) linearly, then there will be 
a linear change of photo-strain, that is new natural state 
of the cantilever, through the thickness. If the geometric 
bend follows this exactly, then there is no effective in- 



ternal stresses in the material. This condition fixes the 
curvature as d t /R = A/a 2 . We thus see that this defi- 
nition of very strong irradiation involves the thickness w 
in relation to adt ■ It is a condition that penetration is so 
deep that light emerges from the other side of a heavily 
dye-loaded sample since there is no region of exponen- 
tial decay of X{x). Moreover, because of the deep pen- 
etration, the curvature now diminishes with increasing 
intensity (like 1 /a 2 ), rather than increasing. The time 
taken for the profile to reach the highly bleached form 
above can be long, a problem we return to in modeling 
the non- linear dynamics of intense absorption [l9| . 

Quantitative measurements of reduced curvature w / R 
with I Q are required to probe this complex dependence of 
curvature on thickness and incident intensity. Particular 
care must be taken to irradiate long enough for equi- 
librium to be reached, and possibly curvature to reduce 
from high intermediate values. 

V. STRAIN DISTRIBUTIONS 

In the linear case [l3[ of bend induced by exponentially 
decaying optical intensity, two neutral surfaces arise, that 
is surfaces of zero stress where the geometric strains aris- 
ing from curvature happen to match the local photo- 
strain: x n /R + K + n c (x n ) =0. A classical cantilever 
bent by imposed terminal torques has a single neutral 
surface at its mid-point x n — w/2, so that stresses that 
are equal and opposite about the x = w/2 sum to zero 
to give no net force. In the current case of cantilevers 
bending because of strains generated internally by light, 
rather than by external imposition of torques, the addi- 
tional constraint of no net torque leads to a more complex 
distribution of stresses which gives rise to more than one 
neutral surface. 

The maximum values for both the curvature 1/R and 
contraction K arc intimately related to the positions of 
the neutral surfaces. Differentiating eqns (fTfj)) and (fT7|) 
with respect to the cantilever thickness w and solving be- 
tween the resulting equations, we obtain the relationship: 

<m _ zif_dK_ 

dw w dw ' 

thus both 1/R and K are maximized for the same can- 
tilever thickness. Furthermore, returning dR/dw — 
dK/dw = to either of the differentiated equations, one 
finds that the thickness at which this happens, w*, sat- 
isfies w* /R + K + n c (w*) — 0, i.e. the straight line of 
geometrically induced strain intersects the photo-strain 
(~ n c (x)) curve at the back surface of the cantilever. 
The back surface is a neutral surface when the curva- 
ture is maximized. The curve n c {x) is always a strictly 
decreasing function of x; however there is a fundamen- 
tal difference between the curves in the (3 — case for 
a < 0.5 and those for a > 0.5 - the latter n c (x) curves 
have a point of inflection (n"(x) — 0). The straight line of 
geometrically generated strain imposed from curvature, 
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FIG. 7: Neutral surfaces for fixed incident light intensity a = 
5, a//3 — 50 as cantilever thickness w changes. At the w 
maximizing curvature, a third neutral surface appears from 
the rear face; at greater w a neutral surface is lost at the front 
face. 



— (x/R + K), can intersect the n c (x) curve at most twice 
if there is no point of inflection, and at most three times 
if there is a single point of inflection. It is not possible 
for the strain to satisfy the constraints of vanishing force 
and torque, that is satisfy eqns (fT5)) . and have a neutral 
surface at the back surface of the beam with only two 
neutral surfaces; thus there is no maximum for the cur- 
vature unless the underlying cis curve has an inflection 
point. For /3 ^ the a value that divides the curves into 
those with and those without an inflection point is no 
longer exactly a = 0.5 but, as can be seen in Fig [31 for 
small (3 the division is still at an a close to 0.5. 

To illustrate the significance of inflections, two sample 
curvature strains of purely geometric origin are superim- 
posed in Fig [3] On the a = 10, j3 = 0.2 curve is also 
plotted the straight line —(x/R + K) for a cantilever of 
thickness w = 12.5dt, a thickness that is before the in- 
flection point in the n c (x) curve. As can be seen, the 
straight line intersects the underlying n c (x) curve only 
twice in satisfying eqns (|15p and locates only two neutral 
surfaces. On the a = 5, f3 = 0.1 curve is also plotted 

— (x/R + K) for the critical thickness w — w* — 9.313<it 
where the third neutral surface first appears. One sees 
that this line intersects the underlying n c (x) curve three 
times, internally twice with the third intersection (neu- 
tral surface) coinciding with the back surface. This line 
is that of maximal possible slope, dR/dw = 0. With 
further increasing thickness, d/R decreases. Eventually 
a neutral surface migrates to the front face and is lost. 
The cantilever continues to have only 2 neutral surfaces 
thereafter. Fig. [7] shows how the neutral surfaces change 
with increasing thickness at fixed illumination a = 5, 
a/ '/3 = 50. Features are slightly shifted to smaller w in 
comparison with the (3 = case of [l6| . With two regions 
of compression and two of elongation, one expects subtle 
behavior when considering compression-induced director 
rotation [2~i| in such cantilevers, to which we return else- 
where. 



VI. TEMPERATURE EFFECTS 

We have neglected the effect of heat generated by ab- 
sorption of light. Gradients of optical intensity in the 
cantilever might be expected to generate gradients of 
temperature and thus of thermal contraction, leading 
to thermal bend. Experiments where polydomain elas- 
tomers bend in the direction of the polarization of light 
[Il| show directly that optical effects dominate (see the 
analysis of polydomains in [Trover thermal component 
of bend. We here quantify the relative size of optical and 
thermal effects. Let the temperature distribution in the 
beam be 6(x,t), and take the origin of the temperature 
scale be such that the ambient temperature outside the 
cantilever is zero, 9 = 0. The temperature distribution 
satisfies a continuity equation in which the heat flux con- 
tains the usual thermal gradient term — tffx and the 
divergence of the Poynting flux, I(x)x, is a source term: 



dt dx 2 



-Jo— 

dx 



(24) 



where C is the specific heat of the the cantilever per unit 
volume and k is the thermal conductivity perpendicular 
to the director, that is along the normal to the flat surface 
of the cantilever. The diffusion co-efficient D is given by 
the ratio n/C. Typical values for elastomers are around 
~ 10 _7 m 2 s _1 [25], with some anisotropy in D arising in 
nematic elastomers from anisotropy in the conductivity 
that we ignore here since thermal effects will in any case 
turn out to be small. The time taken for heat to diffuse 
across the thickness of the cantilever is ~ w 2 /D ~ 0.001s 
for a 10/im sample (Tll.[26j. For times significantly longer 
than this, one obtains the steady state solution 9(x). At 
the front and back surfaces there are convective losses 
which are described by Newton's law of cooling, that is 
the heat flux carried away from a surface is 58(0) (the 
temperature outside is 8 = 0), where S is the convective 
heat transfer co-efficient. For free convection in air, 5 « 5 
W iaT 2 K -1 . These convection losses are equal to the 
thermal flux of heat at the respective surfaces: 



8(0)- ^ 
o dx 



x=0 



8(w) 



ndd 
6 dx 







0, 



(25) 



where the signs reflect the direction of the outward sur- 
face normal. The thermal conductivity is k w 0.2 W m _1 
K _1 [27j ■ The solution to eqn (j2"4")l satisfying boundary 
conditions is: 



8(x) 



1 



(l + — 



1 



5w 

K 



l(wy)dy 



Sw . 
— / l(wy)dy 



(26) 



The characteristic temperature 8 = Iq/S is the temper- 
ature that the sample surface would need to attain in 
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order to lose by Newton cooling all the heat equivalent 
to the incident Poynting energy. Both of the integrals 
appearing in this expression have been scaled such that 
their value is bounded by unity. The scale of their contri- 
butions is thus set by their pre- factor Sw/k. This dimcn- 
sionless quantity compares the thickness of the sample w 
with the thermal penetration length set by the boundary 
conditions k/S. 

There are several interesting limits to this equation: 
(i) Using the values for k and 6 given above and 
assuming a cantilever thickness w ~ lOfim we obtain 
— w 2.5 x 10 -4 << 1, and thus the temperature dis- 
tribution is essentially constant throughout the sample, 
that is: 



6{x) 



I -In 



+ 0(6w/k) 



(27) 



A constant temperature through the sample leads to 
contraction along the long axis of the cantilever, but it 
will not induce bending, since this requires differential 
contractions. Therefore in this limit, which is close to 
experimental reality, bending effects in the steady state 
are due to the optical effects discussed previously, rather 
than heating. 

(ii) Heat is produced proportionately to the rate op- 
tical intensity diminishes. There are regimes in which 
intensity T{x) decays linearly with depth into the can- 
tilever, such as the a = 10, /3 = case of 16[ where there 
is linearity throughout if w < 10d t , and here in Fig [2] for 
a = 10, (3 = 0.2 where the range of thicknesses where the 
entire profile is linear is slightly reduced by the influence 
of stimulated back reactions ((3 ^ 0), say w < 8dt- For 
these conditions, that is where dT/dx = const., then heat 
is generated at the same rate through the cantilever. In 
the steady state it diffuses to the surfaces, symmetrically 
if the surfaces are at the same temperature, and hence 
heat does not contribute at all to the bending. One sees 
this since now eqn l|26p becomes: 



6{x) 



I -In 



Sw 
An 



_5_ 

KW 



(*-!)' 



(28) 



a temperature distribution which is indeed symmet- 
ric about the mid-point of the cantilever. Since the 
thermally-induced strains are also symmetric, they will 
not produce bending, although once again they will pro- 
duce an overall contraction. If the cantilever thickness 
w is increased beyond the linear profile interval in e.g. 
Fig [2l then asymmetry in the heat production starts to 
occur, with less heat generated towards the back face. 
However, the magnitude of the asymmetry is reduced by 
diffusion and the thermal contribution to bend sets in 
only slowly with increasing w, see (i) above. 

The neglect of heat is thus justified in two limits, firstly 
the convective heat losses from the boundary are likely 
to result in a uniform temperature distribution through 
the sample for experimentally realistic values of the ther- 
mal conductivity k and the convective heat transfer co- 



efficient S. Further, in the regime of linear (i.e. non- 
exponential) decay of intense beams there is no thermal 
component of bend. 

Thermal effects become more extreme, especially in 
elastomers, if an interior region of the cantilever's tem- 
perature exceeds the nematic-isotropic transition tem- 
perature. At this temperature an extremely large strain 
can develop in elastomers, and possibly in glasses. If 
it occurs in a region symmetrically disposed about the 
cantilever mid plane, it could lead to pronounced con- 
tractions, but still not to bend. If it occurs in an asym- 
metrically disposed region, it could lead to large bends 
- an extreme case we return to in considering specifi- 
cally elastomers and thus at the same time considering 
director rotation. Additionally one would there consider 
the displacement of neutral planes due to volumes of 
the cantilever suffering large contractions of their nat- 
ural lengths. 



VII. CONCLUSIONS 

We have shown that non-linear absorption (that is, 
non-exponential profiles) can be invoked to explain how 




FIG. 8: (a) Scaled intensity at the back surface I w as a func- 
tion of the incident intensity a for various cantilever thick- 
nesses w/dt with j3 = (solid lines), and with a//3 = 50 
(dashed lines) and a//3 — 5 (dotted lines) for two marked 
thicknesses, (b) Intensity plotted in the form of the ab- 
sorbance — ln(2ur) as a function of a with cantilever thick- 
nesses and p values as in the upper figure. 
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bending can arise in cantilevers where, within Beer's law, 
one would otherwise expect no response. At high enough 
incident light intensities, there can be photo-bleaching 
and thereby penetration of radiation and thus elastic 
response even in cantilevers so heavily loaded with dye 
that the Beer penetration depth of the linear regime is 
insignificant compared with the cantilever's thickness. 
The non-linear response in the high dye-loaded limit is 
possibly of the greatest experimental relevance. We have 
explored the roles of optically-stimulated de-excitation 
of dyes and of optically-generated heat in mechanical 
processes. Stimulated decay can be very important, but 
the role of heat seems minor. 

An experiment [13] to explore the non-linear absorption 
processes we have described measures the intensity 
emergent at the back surface for a sample of fixed width 
w illuminated at the front surface. Altering the incident 
intensity is equivalent to varying the parameter a. 
Taking a fixed depth (the thickness itself) w/dt in Fig [5] 
and [5] and increasing reduced incident power a — Iq / I t is 
to take a slice through these figures to reveal increasing 
relative penetration, that is an increased reduced exit 
intensity T w {a) — I w /Iq, see Fig |5{a). Scaling the 
output intensity by its incident value would produce 
a horizontal line X w — exp(—w/dt) as a function of a 
for simple Beer law attenuation; deviations from this 
line as power increases away from a — are thus signs 
of non-Beer attenuation processes. In thin cantilevers 
w/dt — 0.1, T w is close to unity, and only increases 
slowly as a is increased - for thin cantilevers nearly all of 
the incident flux is transmitted in any case. Conversely 
for thick cantilevers w/d t = 10 we see that there is very 
little transmittance for small a. For larger incident 
intensities, such that a > 5, the transmittance begins to 
increase rapidly with a - a consequence of the increased 
penetration due to non-linear absorption processes. 
For two particular sample thicknesses, |3 / cases 
are shown too. The curves show that a finite induced 
back reaction rate reduces enhanced penetration, as one 
would expect. For the thickest sample, the higher (3 
value is sufficient to effectively eliminate penetration of 
the beam. Deviation from Beer absorbance is perhaps 
better seen from plots of non-linear absorbance itself, 
see Fig. Efb). An experimental example, and its relation 
to their equivalent of our eqn ©, is discussed at length 
by Statman and Janossy pTj . 

Despite host and dye absorptive processes being mixed 
in a non-linear differential equation, an approximate al- 
lowance by division by the host transmittance seems to 
work well: C w = T w /e~ w l dh by substitution in eqn Ijlip 
yields: 



ln(C u 



1 



In 



1 + ^(1 + . 



1 + ^(1 + a) 



-u>/dt.(29) 



We plot in Fig. [9] the corrected exit intensity C w implied 
by ([21]). For very intense beams, a 3> 1, the corrected in- 
tensity tends to 1. The dye is bleached and absorption is 
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FIG. 9: Corrected intensity at the back surface C m as a func- 
tion of the incident intensity a for various cantilever thick- 
nesses w/dt. In each case the upper curve is without host 
absorption, the lower with rather heavy host absorption such 
that w/dh = 0.3. For comparison, lower lighter lines, the 
corresponding uncorrected T w curves are shown. 



predominantly by the host. The true profile is then Beer- 
like and the correction is accurate. At lower intensities 
the correction is also seen to work surprisingly well. The 
upper of each curve for a given reduced thickness w/dt is 
for no host absorption, representing the pure system be- 
ing examined. The corrected curve (lower in each case) is 
close to the ideal curve, even for the very heavy absorp- 
tion w/dh = 0.3 in the illustration. For smaller absorp- 
tions the corrected and the ideal curves are practically 
indistinguishable. Thus the straightforward correction 
method should give a good estimate of the underlying 
non-linear absorptive processes. 

Much of the discussion has been of intensities with ref- 
erence to the ideal Beer intensity, for instance Fig[5fa). 
Beer behavior obtains when intensities are low enough 
that the trans population is little reduced. This also ob- 
tains at short times when as yet little conversion has 
taken place. An output power X w = exp(— w/dt) must 
emerge and has the value of the vertical axis intercept in 
Fig (Ha) (and is also the same as the reduced weak-beam 
emergent power). Thus, for a given power a and thick- 
ness w/dt, the initial emergent power must rise to the 
higher value characteristic of the curve in question. The 
dynamics by which a spatial bleaching pattern is estab- 
lished is non-linear and complicated, see comments below 
eqn 

One can also propose that bleaching can lead to an 
emergent power that diminishes in time from the initial, 
Beer value if the cis species have a higher absorption than 
the trans . A higher cis absorption arises if the incident 
light frequency is shifted to be closer to the cis absorption 
line. One can rearrange eqn ^ by evaluating it at x = w 
and noting that —w/dt = In (2b), where b denotes the 
Beer, t = reduced intensity. Thus denoting by 1Z the 
ratio Zuj/Zb, one has: 



(i + v)/v 



1+0' 



(30) 
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FIG. 10: Intensity at the back surface scaled by the Beer in- 
tensity expected at the back face, 1Z = I w /Ib as a function 
of the incident intensity scaled by the intensity characteristic 
of the as state, that is j3 = Jo /J c for cantilever thicknesses 
w/dt = 4, 1 and 0.5 with j3/a = 50 (solid) and (3/a = 5 
(dashed) pairs of lines for each thickness. The emergent in- 
tensity is lower than the Beer value, in contrast to the rest of 
the paper where bleaching enhances penetration. 



It is evident that the In on the right hand side is al- 
ways negative since TgR = T w < 1 and hence if a < rjf3 
then the left hand side is also negative, that is 7Z < 1 
- the intensity drops below the Beer value rather than 
rising with time and increasing power. The criterion, 
referring back to the ratio of the quantum efficiencies, 
amounts to d c < d tl that is, attenuation due to cis is 
greater than that due to trans . Sample thickness is en- 
coded in 1b = exp(— io/d t ). Solutions to eqn (|30f are 
shown in Fig 1101 Care is needed in interpreting this fig- 
ure: as thickness increases, the absolute amount of light 
passing through the sample is greatly reduced since the 
Beer normalisation to X w in 7Z is becoming exponentially 
small. This reverse effect of photo-bleaching maybe use- 
ful in occasions where at high intensities it is desirable to 
have a smaller fraction of the light transmitted (but as 
noted above, the absolute amount of light transmitted is 
still increasing with incident power). It is possible that 
biological systems depend on this behavior. 
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